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Abstract 


This  Report  discusses  work  on  Grant  AFOSR  86-0119  for  the  period  5/1/90- 
1/31/91.  Previous  work  on  ionization  non-equilibrium  is  summarized,  and 
new  results  pertaining  to  inlet  ionization  in  MPD  thrusters  are  presented. 
Brief  summaries  are  also  given  of  our  work  on  transport  effects.  Hall  effects 
magnetic  layers  and  electrothermal  arejets. 
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This  Final  Report  summarizes  the  results  of  our  work  on  plasma  non¬ 
equilibrium  issues,  under  Grant  AFOSR-86-Ol  19.  Results  to  4/30/90  have  been 
documented  in  our  previous  Annual  Reports,  so  we  will  here  concentrate  on 
the  final  segment  of  the  Grant  covering  the  period  5/1/90-1/31/91. 

The  main  responsibility  for  the  work  on  Radiation  and  Non-equilibrium 
has  been  on  doctoral  candidate  Eric  Sheppard,  who  has  been  directly  supported 
by  the  Grant.  In  addition  the  existence  of  AFOSR's  Grant,  and  the  resultant 
research  activity,  has  attracted  several  other  self-supporting  students,  who 
have  also  contributed  significantly  to  our  progress.  These  include  doctoral 
candidates  J.  Marc  Chanty,  Eli  Nicwood  and  Scott  Miller.  Recently,  the  Air 
Force  Phillips  Lab.  has  entered  into  a  cooperative  program  with  us,  and  a 
doctoral  student,  Mr  Eric  Gaidos,  will  be  starting  experimental  MPD  work  at 
Edwards  AFB  this  Summer.  Also  recently  initiated  is  arejet  modeiing  work  by  a 
new  Master's  student.  Miss  A.  Sakamoto. 

The  object  of  our  research  has  been  the  basic  physics  of  plasma 
thrusters,  particularly  MPD  thrusters.  More  specifically,  within  the  context  of 
the  present  Grant,  we  have  addressed  the  complex  of  issues  which  encompass 
excitation  -  ionization  kinetics  of  MPD  plasmas,  and  their  effects  on  thruster 
flows. 

The  work  on  the  detailed  kinetics  of  electronic  excitation  of  Argon  and 
Hydrogen  was  summarized  in  the  paper  by  E.  Sheppard  and  M.  Martinez- 
Sanchez  (AIAA  90-2608)  which  is  appended  to  this  Report.  Because  of  their 
intimate  coupl  ing,  excitation,  line  radiation  and  ionization  effects  must  be 

considered  together.  The  situation  in  MPD  thrusters  (as  opposed  to  arejets)  is 
far  from  Local  Thermodynamic  Equilibium  (LTE),  and  modeling  has  to  be  based 

on  a  detailed  set  of  atomic  level  population  balances.  We  included  three  basic 
non-equilibrium  drivers. 

(a)  Radiation  escape.  The  majority  of  the  spectrum  is  optically  thin 

(with  the  exception  of  the  central  regions  of  resonant  lines),  and 
there  is  therefore  a  continuous  net  radiative  depopulation  rate  of 

all  excited  states. 

(b)  Pi ffusion.  Due  to  the  prevailing  low  density,  MPD  plasmas  share 

many  to  the  characteristics  of  glow  discharges,  including  the 
prevalence  of  ambipolar  diffusion  as  a  loss  mechanism  for  ion- 
electron  pairs.  This  requires  continuous  excess  ionization  to 
maintain  the  ionic  balance,  and  this  cascades  down  the  excited 

level  populations. 

(c)  Dynamic  Lae.  Any  given  parcel  of  plasma  flowing  through  a 
thruster  experiences  rapidly  changing  background  conditions 
and  many  of  its  internal  degrees  of  freedom  lag  widely  behind 
shifting  equilibrium.  This  is  most  extreme  at  the  channel  inlet, 
where  gas  goes  in  a  few  microseconds  from  room  conditions  to 
over  an  eV  background  energy.  Overall  ionization  equilibrium 
(Saha's  law)  is  clearly  violated.  Other  processes,  such  as 
population/dcpopulation  level  kinetics  remain  close  to  balanced 
in  a  quasi-static  sense. 


Using  multi-level  models  with  the  quasi-static  approximation  for 
internal  balancing,  we  derived  expressions  and  values  for  line  radiation 
density,  non-Boltzmann  population  distributions  and  overall  ionization- 
recombination  rate  coefficients.  It  was  found  that  radiative  loss  can  account 
for  at  most  a  few  percent  of  the  input  power  to  an  MPD  thruster,  in  good 
agreement  with  experimental  (bolometer)  data  from  the  related  field  of  coaxial 
plasma  guns  Ul.  it  was  also  found  that  care  must  be  exercised  when  using  the 
upper  level  populations  of  excited  Argon  ions  for  Te  measurements,  since 
significant  deviations  occur  from  the  normally  assumed  equilibration  with  the 
free  electrons.  This  is  not,  however,  the  case  with  neutral  Argon,  although, 

unfortunately,  its  lines  are  normally  much  weaker. 

2.  Inlet  Ionization  Effects 

The  inlet  relaxation  problem  has  been  one  of  the  most  important  and 
difficult  areas  of  uncertainty  in  plasma  devices,  Non-equilibrium  MHD 
generators  suffer  global  instabilities  due  to  this  problem^2!.  Our  own 
experimental  MPD  thruster^  failed  to  produce  an  azimulhally  uniform 
discharge  when  the  inlet  electric  field  was  too  weak  due  to  a  convergent-type 
design.  Japanese  experiments^  showed  that  Hydrogen  would  ionize  very 
weakly  and  with  large  delays  in  a  thruster  in  which  Argon  showed  no  such 
difficulties.  The  gas  enters  these  devices  with  no  ions  or  free  electrons,  and  a 
rapid  process  of  electron  multiplication  must  take  place,  starting  with  a  few 
back-diffused  or  photoionized  electrons.  As  the  case  of  Hydrogen  vs.  Argon 
shows,  the  details  of  the  kinetics  of  intermediate"radicals"  can  control  the 
whole  process,  in  analogy  with  the  ignition  delay  problems  in  high  speed 
combustors. 

One  type  of  injector  which  is  potentially  important  for  steady-state  MPD 
thrusters  is  a  ceramic  porous  plate  (Fig.  la).  This  is  one  limiting  configuration 

which  can  be  analyzed  as  a  one-dimensional  structure.  In  previous  work,  we 

have  studied  the  dynamics  of  such  an  inlet  layer  assuming  full  ionization.  We 
are  now  extending  this  by  including  kinetic  equations  for  ionization  and 

species  diffusion.  Radiation  and  excitation  kinetics  may  also  play  a  role,  and 
will  be  assessed  later. 

A  more  common  configuration  involves  discrete  injector  holes  (Fig.  IB). 

Here,  the  diffusion  of  electrons  into  the  fresh  gas  is  mainly  across  streamlines, 
rather  than  against  the  flow,  and  the  structure  resembles  a  3-D  diffusion 
flame.  The  study  in  the  three  dimensional  case  is  very  difficult,  and  we  will 
attempt  to  simplify  it  by  using  its  slenderness  to  treat  it  as  a  submerged 
reacting  jet.  No  results  have  yet  been  obtained  for  this  configuratir ,i. 

In  either  case,  the  results  should  indicate  the  length  of  the  inlet 
relaxation  zone,  and  also  the  threshold  conditions  below  which  no  plasma 
ignition  occurs.  Both  of  these  items  are  important  design  information  for 
plasma  devices. 

Significant  results  have  been  obtained  for  the  one-dimensional  case. 
Because  of  the  presence  of  the  insulating  (hence  cicctron-rcpciling)  back 
plate,  the  electron  heat  flux  must  be  small  in  the  ionization  layer,  and  hence  a 
constant  Te  approximation  is  made.  Assuming  also  a  constant  flow  speed 
through  the  layer  (an  assumption  to  be  removed  later),  one  can  integrate  the 
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(a)  Porous-Plug  Injector 

1-0  model  of  Ionizing  front: 


(cold  gas)  plasma 

Questions: 

*  What  processes  control? 

*  How  thick  is  the  front? 

*  Under  what  conditions  is  there  a 
solution?  (ignition  conditions) 


Fig.  1(a).  Schematic  variation 
of  plasma  quantities  near  a  1-D 
(porous  plug)  injector 


(b)  Individual  Injector  Jet 


Fig.  1(b).  The  flow  field 
around  an  individual 
injector  jet. 


Questions: 

*  How  does  backjround  plasma  Interact  with 
the  fresh  cold  gas? 

*  How  long  is  the  “jet"? 

*  How  does  the  “Jet“  interact  with  the 
Lorentz  forces? 
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U  (m/;s) 


concection-diffusion  -  ionization  equations  between  a  wall  condition  relating 
diffusive  ion  flux  to  Bohm  velocity  times  local  density  and  and  a  downstream 
asymptotic  plasma  state.  It  is  found  that  no  attached  ionization  layer  can  exit 
beyond  a  flow  velocity  (which  increases  strongly  with  electron  temperature). 
When  an  attached  ionization  layer  does  exist,  its  thickness  varies  as 

(velocity)®'^  and  decreases  exponentially  with  Te.  ForTe>leVand  reasonable 
conditions,  in  Argon,  the  layer  thickness  is  found  to  be  a  few  mm.  or  less.  The 
near-back  plate  electron  density  is  found  to  be  high  enough  to  support  strong 

electron  thermal  conductivity,  hence  validating  the  Te  =  constant  assumption. 
These  results  are  summarized  in  Fig.  2. 

The  above  model  concentrated  on  the  counteracting  effects  of 
convection  and  diffusion,  and  ignored  the  flow  dynamics.  At  the  opposite 
extreme,  our  earlier  modclsl5]  studied  the  flow  dynamics  in  the  acceleration 
region  while  ignoring  ionization  effects. 


In  reality,  plasma  acceleration  occurs  simultaneously  with  the  Ohmic 
dissipation  which  feeds  ionization,  and  both  effects  should  be  considered 
together.  Retaining  the  constant-temperature  approximation,  ignoring  bulk 
recombination,  lateral  diffusion  and  area  variation  in  the  ionization- 
accerleration  layer,  the  governing  equations  are 


(ne  +  nn)  u  =  G 


m;  G  u  +  ne  k  (Te  +  Tg)  +  nn  k  Tg 

-f  (neu-Da^k)  =  RSn„ne 
dx  dx 

e  Vi  R  S  n„  ne  =  o  (E  -  u  B)2 


=  -  |i0  a  (E  -  u  B) 


(1) 

(2) 

(3) 

(4) 

(5) 


where  ne,  n„  are  the  electron-ion  and  the  neutral  density,  respectively,  m;G  is 
the  known  mass  flux,  F  is  the  (unknown)  momentum  flux,  and  separate 
electron  and  heavy  species  temperatures  (Te,  Tg)  are  assumed,  both  taken  as 
constants. 


In  Eq.  (3),  Da  is  the  ambipolar  diffusion  coefficient  (Da  (ne  +  n„)  =  const.), 
R(Te)  is  a  Saha  factor,  such  that  RS  is  the  ionization  rate  constant  at  the 
prevailing  Te>  This  Te  is  to  be  inferred  from  the  simplified  electron  energy 

balance  in  Eq.  (4).  Since  S(TC)  is  an  exponential  function  ofTe*sma11 
variations  are  to  be  expected,  and  a  mean  value  suffices  for  non-sensitive 
functions  of  Te,  The  one  sensitive  term  is  the  right  hand  side  of  Eq.  (3),  where 
the  whole  source  term  is  to  be  substituted  from  (4).  Finally,  in  Eq.  (5)  the 
conductivity  o  is  assumed  to  be  constant  throughout,  since,  as  will  be  seen, 

electron  diffusion  can  maintain  ionization  fractions  of  over  =  1%  right  down 
to  the  edge  of  the  clcctron-rcpclling  sheath  on  the  injector  plate. 
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Reorganizing  Eqs.  (3)  and  (5),  they  can  be  simultaneously  solved  for  a 

set  of  three  first  derivatives  such  as  and  where  rD-Da^^-is  the 

backwards  diffusive  flux,  when  this  is  done,  it  is  found  that  the  set  of 
differential  equations  has  a  singularity  at 


u  =  u 


* 


(6) 


Local  analysis  then  shows  that  the  singularity  is  of  the  repelling  node 
type,  such  that  no  smooth  crossing  through  the  velocity  u*  is  possible  within 
the  region  described  by  the  given  equations.  Either  the  flow  enters  this 

region  supercriticallv  (u  >  u*)  or,  at  least,  critically  (u  =  u*).  For  Argon  at  a 
few  thousand  degrees  Kelvin,  u*  is  a  few  hundreds  of  m/sec.  How  the  flow 
can  reach  the  entrance  speed  is  related  to  the  constant-temperature 
assumption  introduced,  which  implies  a  steady  supply  of  thermal  energy  by 
conduction  from  regions  downstream.  This  can  provide  the  energy  needed  to 
"thermally  choke"  the  arriving  neutral  gas  flow. 

Thus,  the  condition  u  =  u*  is  a  natural  specification  at  x  =  o,  together 
with  B  =  B0-  Another  one  states  that  the  total  flux  of  ions  at  the  sheath  edge 
must  be  backwards,  and  must  equal  the  sheath-entering  Bohm  flow: 


(n'u'D*frLo“‘nc«VB  ;  VB =V 


(kiy+T.)" 


m; 


(7) 


Finally,  at  the  downstream  end  we  simply  assume  a  length  much  greater 
than  the  ionization-acceleration  distance,  such  that  an  asymptotic  condition  of 
zero  current  and  zero  gradients  be  reached  15K 


This  implies 


du  _  dne  _  dB  _  d£o  _  0  (  x 
dx  dx  dx  dx 

which,  from  Eq.  (5)  requires 


(8) 


The  transverse  electric  field  F.  is  here  arbitrary  (within  limits).  Since 

the  channel  is  thought  of  as  possessing  a  throat  somewhere  downstream, 
different  choices  of  E  will  dictate  different  area  ratios  A*/A0  between  throat 
and  inlet.  The  largest  choice  of  E  is  that  for  which  A*/A0  =  1,  and  smaller  E 
values  correspond  to  A*/A0  <  1- 

For  numerical  integration,  the  behavior  at  large  x  is  important.  An 

asymptotic  analysis  shows  that,  when  the  equations  arc  linerarized  for 

conditions  near  u„,  B^,  etc,  two  characteristic  modes  occur,  of  which  one 
diverges  violently  in  the  downstream  direction,  while  the  other  converges  at  a 
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moderate  rate.  It  is  therefore  necessary  to  integrate  backwards.  The 

procedure  is  roughly  as  follows:  a  set  of  values  ne^,  is  guessed  at.  This  plus 
the  conditions  at  infinity,  gives  sufficient  information  to  calculate 
analytically  the  amplitude  of  the  acceptable  downstream  mode,  and  therefore 
to  "step  out”  of  the  asymptotic  limit  by  some  small  amount.  From  there, 
numerical  integration  is  continued  back  to  the  inlet  condition  B  =  BQ.  At  this 
point  the  conditions  u  =  u*  and  (flux  =  -  Bohm  flux)  are  checked,  and  iteration 
on  B  ,ne  follows  until  they  are  both  met. 

oo  l  oo  * 


Results  for  a  case  with  convergence  ratio  H*/Ho  ~  0.728  and  with  a 
slightly  supercritical  inlet  velocity  are  shown  in  Figs.  3  and  4a-4d.  Several 
dimensionless  parameters  need  to  be  defined  for  interpretation: 


—  ^ref 

ucrv 


Bp 

2|i0  miG 


,  UdV 


0  = 


Tj 

Te  +  T, 


£d 


-g-%  Da  (ne  +  nn) 
2  mi  G2 


(10) 

(11) 

(12) 

(13) 


of  these,  X  measures  the  more  familiar  (current)^/(mass  flow  parameter,  and 
should  be  of  order  1  for  operation  near  "onset”  (A.  =  1  taken  for  our  example). 

0  must  be  small  near  the  inlet,  and  we  use  0  =  0.1  here.  Similarly,  Eq 
measures  the  ratio  of  ambipolar  to  magnetic  diffusivity  (at  the  characteristic 
MPD  density),  and  is  a  small  parameter  (0.1  used  here).  Finally,  p  is  a  measure 
of  the  gas  pressure  divided  by  Magnetic  pressure,  and  the  value  p  =  0.05  has 
been  assumed. 

Fig.  3  shows  the  normalized  magnetic  field  B'  =  B/B0  versus  normalized 
distance 

X  =  O  Ho  Uref  x  (14) 

which  can  be  interpreted  as  the  Magnetic  Reynolds  number  for  a  distance 
equal  to  x.  B'  is  taken  to  decrease  from  1  to  0.9,  corresponding  to  the 
prescribed  inlet  velocity.  The  characteristic  distance  for  B  variation,  namely, 

for  the  initial  current  concentration,  is  seen  to  be  of  the  order  A  x'  =  1-  This  's 
found,  within  a  factor  of  2,  for  all  conditions  locked  at  so  far,  and  is  basically 
the  same  result  obtained  in  Ref.  5  for  frozen  flow. 
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MPD  l-D  INLET 


MPD  l-D  INLET 
RK4S 


Fi.,  4a  B' 
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MPD  I  D  INLET 
RK45 


29  May  91  10:41  30 

Fig  44  B' 
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Fig.  4a  shows  the  dimensionless  velocity  u  =  u/uref  versus  magnetic  field 
B\  Here,  u'  increases  from  the  slightly  supercritical  value  of  0.1  (u*=  0.0707  in 

this  case)  to  u.,  =  0.333,  corresponding  to  =  E/B„,.  Notice  the  nearly  linear 
variation. 


Fig.  4b  gives  the  ionization  fraction  a  vs.  B'.  The  wall  value  a0  =  0.024 
results  from  satisfaction  of  the  boundary  condition  of  Eq.  (7),  and,  as  noted,  is 
sufficiently  high  to  assume  a  nearly  constant  (Coulomb  dominated) 
conductivity  throughout.  The  ionization  fraction  increased  monotonically  on 
approximately  the  same  length  scale  as  u'.  This  is  because  ionization  is  driven 

by  dissipation  (j2/o),  while  u'  is  driven  by  Lorcntz  force  (jB),  both  of  which 
concentrate  in  the  same  region. 

The  normalized  electron  density 
B2 

ne  = - r  ne 

2^0miG2  (15) 

is  shown  in  Fig.  4c.  The  intermediate  maximum  is  characteristic  of  other  cases 
treated  as  well,  and  is  due  to  the  counteracting  effects  of  continuing  ionization 
on  one  hand,  versus  overall  plasma  expansion  on  the  other. 


Finally,  Fig.  4d  shows  the  normalized  net  ion  flux 

*-o(n'u-D*^) 


(16) 


versus  B'.  Notice  the  negative  inlet  (B'  =  1)  value  of  <J),  in  consistence  with  Eq. 
(7) 


So  far,  no  value  of  Te  has  been  assumed.  As  noted,  the  appropriate  Te  is 
dictated  by  Eq.  (4),  which  contains  the  specific  recombination  rate  function 
R(Te).  Work  is  in  progress  now  to  asses  the  effect  of  using  conventional 
models,  such  as  the  Hinnov  -  Hirschberg  equation,  versus  our  more  detailed 
rate  formulation  of  Paper  AFAA  90-2608,  on  the  obtained  electron  temperature. 

3.  Other  Activities 

3.1  Transport  Effects  in  MPD.  Scott  Miller  has  completed  development  of 
a  2-D  code  which  accounts  for  the  full  effect  of  viscosity,  electron  and  heavy 
particles  heal  transfer  and  diffusion,  while,  at  the  same  time  fully  coupling  to 
the  electromagnetic  fields.  This  extends  the  work  reported  in  the  paper  AlAA 
90-2606,  in  which  the  magnetic  field  distribution  was  independently 
prescribed.  Results  for  a  channel  with  2  cm.  intcrclcctrodc  separation,  10  cm 

active  length,  with  nn(Argon)  =  4  g/Scc  and  1  =  35  KA  are  shown  in  Fig.  5.  The 
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Two-Dimensional  Viscous  MPD  Flow 
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first  panel  shows  the  current  lines  to  be  only  slightly  convected/distorted, 
which  is  a  result  of  ignoring  the  Hall  effect.  The  velocity  distribution  shows 
boundary  layers  which  nearly  reach  full  development  at  exit,  and  which 
create  substantial  drag  forces.  Finally,  the  gas  temperature  shows  the  classical 
hypersonic  boundary  layer  overshoots,  peaking  at  about  1  cV  near  the 
channel  exit.  The  lack  of  symmetry  in  the  temperature  and  velocity  maps  is 
due  to  the  effects  of  heat  transport  and  expansion  work  of  the  electron  gas  as  it 
drifts  transversely. 

3.2  The  Hall  Effect  and  its  Consequences 

While  the  solutions  just  discussed  incorporated  many  real  gas  effects, 
they  still  ignored  anisotropic  electron  conductivity.  Eli  Niewood  has  been 
working  on  the  inclusion  of  this  important  effect,  which  presents  great 
numerical  difficulties  due  to  the  steepness  of  the  gradients  it  induces.  By 
gradually  refining  the  formulation,  particularly  at  and  near  the  electrode 
boundaries,  the  latest  results  already  meet  all  the  conservation  tests  to 
reasonable  accuracy.  In  these  calculations,  heat  transport  was  included,  as 
well  as  ambipolar  diffusion,  but  viscosity  was  not.  For  the  same  case  as  that  of 
Fig.  5,  the  results  arc  shown  in  Fig.  6.  The  current  lines  are  now  drastically 
distorted,  especially  near  the  anode  (lower  electrode).  The  Hall  effect  causes 

the  current  lines  to  enter  the  anode  at  an  angle  tan -l(i  fp  = -2- B  is  the  Hall 

en^ 

parameter),  which  becomes  large  at  low  electron  density.  Simultaneously,  the 
♦  — » 

Lorentz  j  x  B  force  pushes  the  plasma  in  a  direction  perpendicular  to  the 
current  vector,  which  near  the  anode  points  mainly  away  from  the  electrode 
surface,  and  creates  a  hydrostatic  pressure  gradient  with  a  deep  minimum  at 
the  surface.  Thus,  these  two  effects  reinforce  each  other,  and  are  eventually 
checked  only  by  diffusive  and  ohmic  damping.  The  result  is  the  creation  of  a 

largely  depleted  thin  layer  where  P  can  reach  values  of  the  order  of  10-20, 

while  P  remains  of  order  1  in  the  bulk.  In  a  uniform  plasma,  P>2-3would 
be  expected  to  create  an  ionization  instability  [61,  but  the  large  gradients  and 
the  immediately  adjacent  wall  serve  to  prevent  its  development  in  this  case 
(although,  in  a  sense,  the  very  formation  of  the  depleted  layer  may  be  thought 
of  as  the  final  stage  of  this  type  of  instability). 

From  a  practical  standpoint,  the  essential  consequence  of  this  effect  is 
the  creation  of  a  large  anode  voltage  drop.  This  is  because  of  the  Hall 
component  of  Ohm's  law: 
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(17) 


O  (E  +  u  x  B  +  —£■  V  ne)  =  j  +  j  x  P 

The  density  gradient  term  can  itself  be  hydrostatically  related  to  the 
Lorenz  cross  force.  At  full  ionization  (which,  as  shown  by  the  second  panel  in 
Fig.  6  does  occur  near  the  anode),  we  have 


kTe^  =  jxB 

dy 


and  then  the  Y  -  component  of  (17)  becomes 


(18) 


Ey  =  u  B  + 

7  o 


(19) 


and,  since  the  X  component  of  the  left  hand  side  of  Eq.  (17)  is  near  zero  at  the 
wall,  jx  =  -  P  jy>  giving 


Ey 


u  B  + 


1  +  2  p2 
a 


Jy 


(20) 


Since  P  is  very  large  in  the  thin  depleted  layer,  integration  of  Eq.  (20) 

gives  a  large  concentrated  voltage  drop  occurring  over  distances  of  the  order 
of  one  mm.  This  is  illustrated  in  Fig.  7,  which  shows  a  transverse  profile  in  the 

vicinity  of  the  current  peak  of  Fig.  6.  Here,  AVa=17  Volts  in  the  first  mm., 

while  the  overall  anode-cathode  voltage  is  only  25  Volts.  Incidentally,  this  AVa 

matches  well  with  that  measured  under  almost  identical  conditions  in  our 

[7,8] 

experiments  on  the  joint  MIT  -  R&D  program 

It  is  appropriate  to  remark  that,  to  our  knowledge,  this  is  the  first 
satisfactory  explanation  ever  reported  for  the  anomalous  anode  drops 
occurring  in  MPD  thrusters  at  high  currents.  Previous  speculation  had 
centered  on  sheath  effects,  perhaps  complicated  by  magnetic  effects.  Our 

results  indicate  clearly  a  Hall  origin  for  these  drops,  which  occur  in  the 
quasineutral  region  of  the  plasma.  This  clarification  may  in  time  open  ways  to 
avoid  such  effects. 

U _ Magnetic  Boundary  Lavers 

The  idealization  in  which  the  Ohmic  component  of  the  electric  field  is 
neglected  when  compared  to  the  back  cmf  is  very  useful  for  understanding 
the  dynamics  of  bulk  plasma,  but  fails  in  regions  where  the  presence  of 
specific  boundary  conditions  impose  the  existence  of  localized  current 
concentrations.  A  prototype  of  this  behavior  is  the  inlet  (and  to  some  extent 
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the  exit)  current  concentrations  in  one-dimensional  accelerators,  treated  in 
Ref.  5.  J.  Marc  Chanty  has  developed  the  theory  of  these  "magnetic  layers"  for 
several  additional  situations: 

(a)  Electrode  layers,  occurring  when  plasma  flows  along  a  flat  wall 
which  transitions  from  insulating  to  conducting. 

(b)  Insulator  layers,  for  the  reverse  transition. 

(c)  Comer  flows,  when  an  electrode  has  a  sharp  convex  comer 
which  causes  the  plasma  to  expand  around  it  to  a  vacuum. 

Cases  (a)  and  (b)  are  tractable  by  means  of  a  similarity  transformation 
closely  related  to  those  familiar  from  viscous  boundary  layer  theory.  This 

allows  the  model  to  include  viscous  and  other  transport  effects  simultaneously 
with  the  magnetic  diffusion  effect.  Exact  solutions  have  been  generated  and 
will  be  presented  at  the  next  Electric  Propulsion  Conference,  Case  (c)  is  much 

more  difficult,  due  to  the  two-dimensionality.  For  small  comer  turning 
angles,  linearization  about  frcc-strcam  is  possible,  and  a  complete  solution  has 

been  obtained  by  the  use  of  Fourier  techniques.  This  shows  the  magnetic 
diffusion  smearing  the  otherwise  sharp  corner  concentration,  and  also 
smoothing  the  current  "fan"  for  distances  of  many  magnetic  interaction 
lengths.  This  solution  will  also  be  presented  at  the  EP  conference.  Work 
continues  for  more  realistic  large  angles,  motivated  by  the  need  for  proper 
modelling  of  the  anode  exit  corner  flows. 

3.4  Arcjct  Modeling 

In  anticipation  of  the  next  AFOSR  Grant,  which  focuses  in  its  first  year 
on  electrothermal  arcjct  performance  prediction,  we  have  initiated  two 

activities  in  this  area: 

(a)  The  MPD  2-D  model  of  Scott  Miller  has  been  expanded  by  addition 
of  extensive  thcrmochemical  and  transport  submodels  for 
hydrogen  (molecular,  atomic  and  ionized),  and  adapted  to 
axisymmctric  geometry.  The  magnetic  field  dynamics  has  been 
replaced  by  an  elliptic  potential  equation  describing  current 
conservation.  We  expect  this  code  to  become  operational  very 
soon,  and  it  will  be  used  for  an  rigorous  analysis  of  the  internal 
arejet  flows  and  a  prediction  of  their  performance.  Extensions  to 
ammonia  and  other  nitrogen  -  hydrogen  mixtures  are 
contemplated. 

(b)  A  new  Research  Assistant,  Miss  Atsuko  Sakamoto,  will  start  a 
parallel  program  in  which  she  will  collect  and  correlate  the 
existing  data  on  arcjccts,  and  will  also  develop  a  lower  level 
(finite  volume)  arcjct  model  to  be  used  to  rapid  performance 
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estimation  and  correlation.  This  model  will  be  tuned  by 
comparison  to  both,  the  data,  and  the  results  of  the  rigorous  2-D 
calculations. 


3.5  MPD  Experiments  at  Edwards  AFB 

Doctoral  candidate  Eric  Gaidos  will  be  initiating  this  Summer  a  joint 
experimental  program  with  the  Phillips  Lab.  The  goal  is  the  investigation  of 
electrothermal  stability  of  the  plasma  near  full  ionization,  in  an  attempt  to 
verify  the  theory  of  Refs.  9  and  10.  Funding  is  provided  in  part  by  the  NSF 
(Doctoral  Fellowship)  and  in  part  by  the  Phillips  Laboratory. 
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Appendix 


Publications  during  the  Reporting  Period 

Included  arc  the  papers 

•  E.H.  Niewood  and  M.  Martinez-Sanchez.  "Quasi  One-Dimensional 
Numerical  Simulation  of  MPD  Thrusters".  Accepted  for  publication  at 
the  AIAA  J.  of  Prop,  and  Power. 

•  E.  Sheppard  and  M.  Martinez-Sanchez.  "Non-Equilibrium  Ionization  in 
Plasma  Accelerators".  AIAA  90-2608,  july  1990. 

•  J.  Preble  and  M.  Martinez-Sanchez  "Onset  in  MPD  Thrusters:  A  Model  of 
an  Electrothermal  Instbility".  Submitted  to  the  AIAA  J  of  Prop,  and 
Power. 

•  S.  Miller  and  M.  Martinez-Sanchez.  "Viscous  and  Diffusive  Effects  in 
MPD  Flows."  AIAA  90-2606  July  1990. 
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